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The elastic scattering amplitudes of indistinguishable, bosonic, strongly-polar molecules possess 
universal properties at the coldest temperatures due to wave propagation in the long-range dipole- 
dipole field. Universal scattering cross sections and anisotropic threshold angular distributions, 
independent of molecular species, result from careful tuning of the dipole moment with an applied 
electric field. Three distinct families of threshold resonances also occur for specific field strengths, 
and can be both qualitatively and quantitatively predicted using elementary adiabatic and semi- 
classical techniques. The temperatures and densities of heteronuclear molecular gases required to 
observe these univeral characteristics are predicted. PACS numbers: 34.50.Cx, 31.15.ap, 33.15.-e, 
34.20.-b 



Observations of collective phenemena in BEC's of 
magnetically dipolar atomic Cr :l| have spurred at- 
tempts to cool and trap heteronuclear molecules with 
large electric dipole moments. There has been much 
recent progress; for example, 3 x 10 3 cm~ 3 OH [2] 
molecules have been cooled through Stark deceleration 
and trapped at 30 mK; 10 9 cm -3 vibrationally excited 
RbCs |3j molecules, formed via photoassociation, have 
been trapped at 250 [iK\ and 10 12 cm~ 3 (fermionic) 
KRb [4] molecules in ground ro-vibrational states have 
been trapped at 350 nK . A wide variety of techniques 
are now being developed and applied to produce quan- 
tum degenerate gases of such strongly dipolar molecules, 
in search of novel collective phenomena and for various 
applications Q. 

One might naturally expect that as dilute polar gases 
are cooled toward the absolute zero of temperature, a 
regime should emerge in which only the longest ranged 
forces - dipole-dipole interactions - are relevant; in which 
case these different polar gases would share a univer- 
sal equation of state. Myriad theoretical models of ul- 
tracold dipolar gases, most based on mean-field theory, 
assume purely dipolar interactions between molecules, 
modified only by a short-ranged species-dependent con- 
tact potential. [6|. However, recent detailed calculations 
of ultracold molecular collisions show strong species- 
sensitive resonant characteristics of elastic scattering 
cross sections, and display universality only in the high- 
temperature regime of semi-classical collision dynam- 
ics 0. 

One purpose of this Letter is to demonstrate several 
universal characteristics of dipole-dipole collisions that 
will emerge at the lowest temperatures where collective 
phenomena associated with quantum degeneracy are ex- 
pected [§]. Among these universal properties are: 1) 
an absolute minimum elastic scattering cross-section for 
zero-energy collisions of all bosonic dipoles; 2) a strong 
correlation between threshold angular distributions and 
the magnitude of threshold elastic scattering cross sec- 
tions; 3) several sequences of threshold resonant states 
that can be tuned with an applied electric field; and 4) 
a set of effective barrier heights which retard threshold 



scattering in all but one elastic channel. These features 
are independent of short-range interactions between the 
molecular collision partners. 

The near-threshold resonances mentioned above have 
appeared in earlier calculations of RbCs and SrO scat- 
tering [HI , and have also recently appeared in calcu- 
lations of two molecules in a trap [llj . A second purpose 
of this Letter is to provide a complete classification and 
interpretation of these resonances. 

Many molecules possess polar charge distributions in 
their lowest electronic energy states. However, in the 
absence of external fields, as a consequence of parity 
conservation, even strongly "polar" molecules have no 
net dipole moment. Strong anisotropic molecular in- 
teractions therefore require an external electric field, £, 
coupling opposite parity states, and producing a field- 
dependent dipole moment ft(£) = — (dE/d£)£ for all 
molecules 1 101. The magnitude and sign of the "in- 
duced" dipole moment, //, depends on details of the E(£) 
Stark shift of the state of interest for a given molecule. 
Tuning the dipole moment (and therefore the strength 
of the dipole-dipole interaction) with external field per- 
mits variation of the temperature scale at which effects 
of quantum degeneracy can be expected. 

At intermolecular separations much greater than the 
molecular size, molecular interactions reduce to the fa- 
miliar dipole-dipole form 
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V(r) = -o \pi • M2 - 3(f • /2i)(r • p 2 )\ 



(1) 



which, for a pair of identical molecules in the same field- 
aligned state, yields an equation of relative motion 

t 2 



V 2 + ^ (l - 3(f • zf) 



2M 



t/j(f) = E Icl yj(r) (2) 



where r is the intermolecular displacement, z is the field 
axis and M is the reduced mass, fj, assumes its maximal 
value at large fields, but reduces to zero as £ — ► 0. 

Equation © may be written in dimensionless form in- 
dependent of M and /i if all distances and energies are 
measured in dipole units 



D = i?Mjh\ E D = &7(MV) 



(3) 
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The emergence of these field-dependent length and en- 
ergy scales, previously noted by many authors [12], is 
striking when one recognizes that for maximal /i, D ~ 
10 2 — 10 6 times larger than typical molecular length 
scales, while Ed ~ 10 -4 — 10 -12 times smaller than typ- 
ical rotational level splittings. For the examples cited 
in our introduction, we note that D ~ 30 a.u. and 
Ed ~ 10~ 2 K for both triplet KRb and atomic Cr; D is 
of order 10 4 a.u. for OH, singlet KRb, and ground state 
RbCs, while Ed ~ 1 — 100 nK. For a stronger alkali- 
halide dipole such as LiF, D ~ 10 6 a.u. and Ed ~ 0.1 
nK. Note that the recent singlet (fermionic) KRb experi- 
ment is within a factor of 4 in temperature of the intrinsic 
dipole energy scale 0] , and the gas density is comparable 
to l/D 3 . Note also that the energy scales (and relevant 
temperatures) cited here can be increased by decreasing 
the applied electric field. 

There are two restrictions on the domain of validity of 
Eq. ([2]). The dipole-dipole interaction accurately rep- 
resents the intermolecular potential only where higher 
multipoles and exchange terms are negligible. Within an 
exchange radius ~ 20 — 30 Bohr, intermolecular interac- 
tions are complex and species-dependent. There is also 
a range of separations at which the fields created by the 
dipoles dominate the external field, so that effects of colli- 
sional depolarization become important for r ~ [d/S ) 1 ^ 3 , 
typically ~ 100 Bohr in moderate fields for strongly polar 
molecules. Key to the present study is that both of these 
length scales are orders of magnitude smaller than D for 
many systems. 

Equation [2] is usually analyzed by partial wave expan- 
sion [131 ]. In contrast, we adopt here Born-Oppenheimer 
adiabatic states, whose eigenenergies relevant to bosonic 
molecules are displayed in dipolar units in Fig. [T] The 
potential curves and coupling matrix elements of this rep- 
resentation are universal (i.e. independent of collision 
partner, reduced mass, dipole moment, and electric field 
strength) for all r at which the dipole-dipole interaction 
is dominant. Fig. [1] shows that an infinity of channels 
(corresponding to the infinity of partial waves) produce 
a single attractive potential, Vo,0; and an infinite number 
of potentials with centrifugal barriers. 

Due to cylindrical symmetry about the field axis, the 
adiabatic states couple only to other states of the same 
m, which is strictly conserved. The adiabatic curves 
approach familiar centrifugal potentials at large r, per- 
mitting them to be labeled by their large-r angular mo- 
mentum quantum number, £, as well. Accordingly, the 
curves in Fig. [TJ and the resonant states they support, 
are uniquely designated by (I, |m|). Note that V2,i and 
^2,0 present universal barriers of heights ~ 14 Ed and 
~ 140 Ed, respectively. Since m is conserved, the V^i 
channels arc not coupled to the Vo.o channel, while V2,o 
and the ground state channel are strongly coupled. Ac- 
cordingly, at energies much smaller than ~ 140 Ed, scat- 
tering is dominated by the two lowest adiabatic poten- 
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FIG. 1: (Color online) Adiabatic potential curves for bosonic 
polar molecules. The curves are labeled by their asymptotic 
angular momentum £ and by the magnitude of the conserved 
magnetic quantum number m. 

tials, absent the resonances discussed below. 

The most salient feature of Fig. (1) is that the bar- 
rier heights are sufficiently high that tunneling is strongly 
suppressed in all channels near threshold. Accordingly, 
molecules with energies << 14 Ed can approach one 
another to distances less than 0.2D only through the at- 
tractive Vo ; o potential well. While strong non-adiabatic 
coupling exists both inside and outside the potential bar- 
riers, the two regions are effectively decoupled by tunnel- 
ing suppression. It is this feature that leads to universal 
scattering characteristics that are independent of short- 
range physics. This independence is illustrated here by 
applying a variable hard-sphere boundary condition at a 
small radius, r . 

We have solved the coupled-channel Schrodinger equa- 
tion in the adiabatic representation numerically [16j |. 
propagating from ro to asymptotic distances r ~ 5001?. 
While non-adiabatic coupling cannot be neglected, the 
adiabatic basis converges far more rapidly than a partial 
wave expansion near threshold [Hj]. 

The contribution of various channels to the total elas- 
tic cross-section for bosonic scattering (averaged over the 
electric-field direction), is shown in Fig. [2 for ro = 0.0111 
dipole units. The uncoupled Vb,o and V2.1 channels domi- 
nate scattering below ~ 14Ed, as indicated above. Note, 
however, that detailed convergence of the cross-section 
(even at energies less than a nano-Kelvin) requires sev- 
eral adiabatic channels. Both the magnitude of the 
threshold cross section and the degree of convergence 
vary with the hard-sphere radius, as they are sensitive 
to resonance formation in the inner wells of the potential 
curves. This is illustrated in Fig. [3J where the elastic 
cross section (averaged over the field direction) is plotted 
versus hard-sphere radius at the near-threshold energy, 
E = 10~ 3 Ed- Note that the threshold cross-section 
varies by 4-orders of magnitude, depending on the spe- 
cific details of the short-range scattering. 

Figure [3] also illustrates the partial wave content of the 
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FIG. 2: (Color online) Channel contributions to the threshold 
cross section for rO = 0.0111 dipole units. The (0,0) and de- 
generate (2, 1) channels dominate scattering below ~ 14 Ed- 



scattering cross sections. The near-threshold cross sec- 
tion varies quasi-periodically with r due to successive 
"Ramsauer minima" in the attractive adiabatic channel 
(i.e. varying tq alters the number of bound states in the 
ground state potential curve, and correspondingly sweeps 
the s-wave scattering length through complete cycles). 
Despite this sensitivity, the cross section receives a contri- 
bution of 2D 2 independent of Tq from all asymptotically 
repulsive channels. Accordingly, at a minimum of the to- 
tal cross section, the (0, 0) channel is effectively "turned 
off" , and the long-range part of the intermolecular inter- 
action dominates. Absent the resonances, the threshold 
cross section is accordingly governed purely by the long- 
range physics of the dipole-dipole interaction [3]. 




FIG. 3: (Color online) Near-threshold cross section as a func- 
tion of the cut-off radius. Positions of the (2, 0) Feshbach- 
shape resonances are marked by Xs and labeled as Rl, R2, 
R3, Rl 



The resonances shown in Fig. 3 can be predicted by 
semiclassical quantization in accurate fits to the potential 
15| . The positions of the successive Ramsauer 



curves 



\n) = 



1 



4.37 + 3.46?i + 0.675n 2 



(d.u.). 



(4) 



Similarly, the semiclassical description predicts positions 
of resonances in the V2.0 potential that are marked as Rl, 
R2, R3 etc. in Fig. [31 and given by 



2,0 



(n) = 



1 



10. 



II.8I71 + 0.685n 2 



(d.u.). (5) 



If nonadiabatic coupling were negligible, the (2,0) reso- 
nances would have a near-zero (tunneling) width when 
approaching threshold because of the large dipole- 
modified centrifugal barrier in Fig. 1. To the extent 
that tunneling is negligible, these are effectively discrete 
states coupled to the ground state continuum. Their 
actual width is due to decay into the lowest adiabatic 
channel through short-range nonadiabatic coupling. Ac- 
cordingly, the (2, 0) resonances appear as Feshbach/Fano 
resonances in the total elastic cross section. Finally, there 
is also a series of extremely narrow tunneling resonances 
in the V2.1 potentials, whose positions are given by 



2,1 



(n) = 



3.975 + 7.973n + 0.675n 2 



(d.u.). (6) 



but these have little relevance to threshold elastic scat- 
tering, as the |m| = 1 channels are not coupled to the 
ground state channel. 

The dual Feshbach/shape nature of the (2, 0) resonance 
leads to strong correlation between the anisotropy of the 
differential cross section and the magnitude of the to- 
tal cross section. One might expect strongly anisotropic 
scattering near a (2, 0) resonance. However, due to the 
small value of the tunneling amplitude, the inner well of 
the (2, 0) potential is only accessed through coupling to 
the (0, 0) channel. As illustrated in Fig. 4, anisotropy 
results only when the (0, 0) channel is "turned off" (ei- 
ther near a Ramsauer minimum or - due to interference 
- to the side of a Feshback/shape resonance), and stems 
from long-range non-adiabatic coupling outside the (2, 0) 
potential barrier, as well as from the uncoupled (2, 2) 
and (2,1) channels. As a result, whenever the scatter- 
ing cross section peaks, the scattering becomes isotropic, 
even at the peak of the (2, 0) resonance; while near the 
minimum of the total cross section, a strongly anistropic 
distribution emerges. This view augments the interpre- 
tation of similar resonances in atomic collisions in very 



strong fields [17 1 



peaks in the (0, 0) potential are given (to 3 significant 



While the existence of threshold resonances depends on 
ro and will, accordingly, vary from molecule to molecule, 
the resonances can be observed (or eliminated!) for any 
given molecular species by varying the applied field [1 01 ] - 
We described above a dimensionless threshold cross sec- 
tion that varies with a dimensionless cut-off o~a[ro(d.u.)]. 
To the extent that short-range physics is insensitve to the 
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(a) ro = 0.014013 d.u. (b) r = 0.014549 d.u. (c) r„ = 0.014576 d.u. (d) r = 0.014579 d.u. 




FIG. 4: (Color online) Cross section dependence on the angle between the incident direction and the field axis, for values of the 
cut off radius in the vicinity of the Feshbach/shape resonance. When the total cross section is large, the scattering is almost 
isotropic, even at the "d-wave" resonance, (c). 




FIG. 5: Electric field dependence of the near-threshold LiF 
elastic scattering cross section, in both (a) absolute and (b) 
dipole units. The minima of the cross section (= 2D 2 ) in the 
upper figure rise rapidly with the field as a result of rapidly 
growing dipole length. 



external field, the observed elastic scattering cross section 
will scale with the field £ as a — D(£) 2 ao[ro(a.u.) / D(£)], 
where D(£) is given by Eq. [3] This field dependence is 
displayed in Fig. [5] for the LiF molecule, with ro cho- 
sen arbitrarily at 80 a.u., and assuming a quadratic field 
variation of the LiF ground state. The actual, physical 
value of tq can be determined empirically by observation 
of the first peak at the smallest value of the applied elec- 
tric field. This, in turn, lends predictive value to Eq. (J4|, 
which, with knowledge of \i{£ ), predicts the field values 
for all subsequent (0, 0) resonances. The cross-section is 
displayed in both conventional and in dipole units. Also 
shown in the figure is the minimal value of 2D 2 and its 
variation with field strength. 

The adiabatic potential curves serve to identify E < 
lAEjj as the universal energy range of ultracold molec- 
ular collisions, and a ~ 2D 2 as the minimum elas- 
tic cross section for threshold collisions between polar- 
ized molecules. Far greater cross sections approaching 



the unitarity limit ~ 4ttD 2 /(E/E d ) are achievable with 
small variations of the applied field: note the micron 2 
scale in Fig. [5ja) . Observation of the universal proper- 
ties predicted here requires a stongly dipolar gas whose 
temperature is ~ Ejj/k and whose density is less than 
1/D 3 such that two-body collisions dominate. Recent 
experimental efforts are approaching this regime. 

The authors are grateful to John Bohn for helpful con- 
versations and a critical reading of the manuscript. 
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